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Abstract 

Social interactions influence our thoughts, opinions and actions. In this paper, social interactions are studied within a group of 
individuals composed of influential social leaders and followers. Each person is assumed to maintain a social state, which can 
be an emotional state or an opinion. Followers update their social states based on the states of local neighbors, while social 
leaders maintain a constant desired state. Social interactions are modeled as a general directed graph where each directed 
edge represents an influence from one person to another. Motivated by the non-local property of fractional-order systems, the 
social response of individuals in the network are modeled by fractional-order dynamics whose states depend on influences from 
local neighbors and past experiences. A decentralized influence method is then developed to maintain existing social influence 
between individuals (i.e., without isolating peers in the group) and to influence the social group to a common desired state 
(i.e., within a convex hull spanned by social leaders). Mittag-Leffler stability methods are used to prove asymptotic stability 
of the networked fractional-order system. 
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1 Introduction 

Social interactions influence our emotions, opinions, and 
behaviors. Technological advances in social media pro- 
vide more rapid, convenient, and widespread commu- 
nication among individuals, which leads to a more dy- 
namic interaction and influence. For example, recent ri- 
ots [1] and ultimately revolution [2] , have been facilitated 
through social media technologies such as Facebook, 
Twitter, YouTube, and BlackBerry Messaging (BBM). 
Marketing agencies also have begun to take advantage of 
influence due to social media, especially through the in- 
ternet. The company Razorfish, for example, works with 
online peer influencers to transform them into brand ad- 
vocates through the execution of Social Influence Mar- 
keting (SIM) Strategy, which aims to influence market- 
ing primarily through online, small groups, peer pres- 
sure, reciprocity or flattery [3], 

Various dynamic models have been developed to study 
the individual's social behavior, such as the efforts to 
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model the emotional response of different individuals 
[4-6]. In [4], the time- variation of emotions between 
individuals involved in a romantic relationship is de- 
scribed by a dynamic model of love, and in [5] a set of 
differential equations are developed to model the indi- 
vidual's happiness in response to exogenous influences. 
Fractional-order differential equations are a generaliza- 
tion of integer-order differential equations, and they ex- 
hibit a non-local integration property where the next 
state of a system not only depends upon its current state 
but also upon its historical states starting from the initial 
time [7]. Motivated by this property, many researchers 
have explored the use of fractional-order systems as a 
model for various phenomena in natural and engineered 
systems. For instance, the works in [4] and [5] were re- 
visited in [8] and [9] , where the models of love and hap- 
piness were generalized to fractional-order dynamics by 
taking into account the fact that a person's emotional 
response is influenced by past experiences and memo- 
ries. However, the models developed in [4,5,8,9] only fo- 
cus on an individual's emotional response, without con- 
sidering the interaction with social peers where rapid 
and widespread influences from social peers can prevail. 
Other results, such as [10, 11] and the reference therein, 
studied the interaction of social peers using an opinion 
dynamics model, and derived conditions under which 
consensus can be reached. However, agents in [10, 11] 
only update their opinions by averaging the neighbor- 
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ing agent opinions, without taking into account the in- 
fluence of agents' past experience and memory on their 
decision making. 

When making a decision or forming an opinion, indi- 
viduals tend to communicate with parents, friends, or 
colleagues and take advice from social peers. Social con- 
nections such as friendship, kinship, and other relation- 
ships can influence the decisions they make. Some indi- 
viduals (e.g., parents, teachers, mentors, and celebrities) 
may exhibit more powerful influences in others' decision 
making, and the underlying social network enables the 
influence to pass from influential individuals to receptive 
individuals. Containment control is a particular class 
of consensus problems (see [12, 13] for a comprehensive 
literature review for consensus problems), in which fol- 
lower agents are under the influence of leaders through 
local information exchange in a leader- follower network. 
In results such as [14-17], distributed containment con- 
trol algorithms are developed for agents with integer- 
order dynamics where the group of followers is driven 
to a convex hull spanned by multiple leaders' states un- 
der an undirected, directed or switching topology. This 
paper examines how such methods can be leveraged to 
manipulate a social network. This work specifically aims 
to investigate how peer pressure from social leaders af- 
fects consensus beliefs (e.g., opinions, emotional states, 
purchasing decisions, political affiliation, etc.) within a 
social network, and how an interaction algorithm can be 
developed such that the group social behavior can be 
driven to a desired end (i.e., a convex hull spanned by 
the leaders' states). 

By modeling human emotional response as a fractional- 
order system, the influence of a person's emotions within 
a social network is studied, and emotion synchronization 
for a group of individuals is achieved in our recent prelim- 
inary work [18, 19]. However, the emotion synchroniza- 
tion behavior in [18] only considers an undirected net- 
work structure: the one-sided influence of social leaders 
is not considered. This work aims to investigate how the 
social beliefs (e.g., emotional response, opinions, etc.) of 
a group of individuals evolve under the influence of social 
leaders. Similar to [18], the social group is modeled as a 
networked fractional-order system, where the social re- 
sponse of each individual is described by fractional-order 
dynamics whose states depend on influences from social 
peers, as well as past experiences. Since social leaders 
are considered, the undirected network topology in [18] 
is extended to a directed graph, where the directed edges 
indicate the influence capability between two individuals 
(e.g., the leaders can influence the followers' state, but 
not vice versa) . The goal in this work is to develop a de- 
centralized influence algorithm where individuals within 
a social group update their beliefs by considering beliefs 
from social peers and the social group achieves a desired 
common belief (i.e., the social state of the group con- 
verges to a convex hull spanned by social leaders) . Since 
an individual generally only considers others' beliefs as 



reasonable if their beliefs differ by less than a threshold, 
social difference is introduced to measure the closeness 
of the beliefs between individuals. In contrast to the con- 
stant weights considered in [14 16], the social difference 
is a time-varying weight which depends on individuals' 
current states. Moreover, instead of assuming network 
connectivity (i.e., there always exists a path of influence 
between any two agents) such as in [14-16], one main 
challenge here is to influence the followers' social states 
to a desired end by maintaining consistent interaction 
among social peers and influential leaders (i.e., individ- 
uals can always be influenced by social peers, instead 
of being isolated from the social group) within a time- 
varying graph. When modeled as a networked fractional- 
order system, the development of a containment algo- 
rithm can be more challenging compared to the integer- 
order dynamics in [14-17], which can be considered as a 
particular case of generalized fractional-order dynamics. 
The first apparent result that investigated the coordi- 
nation of networked fractional systems is [20]. However, 
only linear time-invariant systems are considered in [20], 
where the interaction between agents is modeled as a link 
with a constant weight. Due to the time- varying weights 
considered here, previous stability analysis tools such 
as examining the eigenvalues of linear time-invariant 
fractional-order systems (cf. [9,20,21]) are not applicable 
to the time-varying networked fractional-order system 
in this work. To address these challenges, a decentral- 
ized influence function is developed to achieve contain- 
ment control for the networked fractional-order systems 
while preserving continued social interaction among in- 
dividuals. Asymptotic convergence of the social states 
to the convex hull spanned by leaders' states in the so- 
cial network is then analyzed via LaSalles's invariance 
theorem [22] , convex properties [23] and a Mittag-Leffler 
stability [24] approach. 

2 Preliminaries 

Consider a Fractional Order System (FOS) 

to V?x(t)=f(t,x) (1) 

with initial conditional x (to), where t 0 X>" denotes the 
fractional derivative operator with order a G (0, 1] on a 
time interval [to^L & n d f (t,x) is piecewise continuous 
in t and locally Lipschitz in x. Similar to the exponen- 
tial function used in solutions of integer-order differen- 
tial equations, the Mittag-Leffler (M-L) function given 

by E ai/ 3 (z) = E £Lo r(fca+/3) > where a, /? > 0 and z e 



2 The initial condition x (to) is defined as a linear combi- 
nation of internal states z^ (to), k — 1, . . . , J, where Zk (to) 
contains all historical information of the system for t < to 
based on the work in [25]. The infinite state model approach 
to resolve the initialization in [25] is also used in the subse- 
quent simulation section. 



2 



C, is frequently used in solutions of fractional-order sys- 
tems. Particularly, when a = f3 = 1, E\,\ (z) = e x is an 
exponential function. Stability of the solutions to (1) are 
defined by the M-L function as follows [24] . 

Definition 1 [24] (Mittag-Leffler Stability) The solu- 
tion of (1) is said to be Mittag-Leffler stable if \\x (t)\\ < 

{m [x (to)] £a,i (—A (t — to) a )} b , where to is the initial 
time, a G (0,1), b > 0, A > 0, m (0) = 0, m(x) > 0, 
m (x) is locally Lipschitz, and E at i is defined in E a ^ 
with [3 = 1. 

Lyapunov's direct method is extended to fractional- 
order systems in the following Lemma to determine 
Mittag-Leffler stability for the solutions of (1) in [24]. 

Lemma 1 [24] Let x = 0 be an equilibrium point for 
the fractional order system (1), and D C K" be a domain 
containing the origin. Let V (t, x) : (0, oo] x D — > TSL be a 
continuously differ entiable function and locally Lipschitz 
with respect to x such that 

<7i(||x||)<K (*,*)< a 2 (\\x\\), 
0 D^V(t,x)<-a 3 (\\x\\), 

where x € H>, f3 S (0, 1) , and <Ji (i = 1, 2, 3) are class K, 
function^]. Then x = 0 is Mittag-Leffler stable, which 
implies that the equilibrium point of (1) is asymptotically 
stable. 

Definition 2 [23] For a set of points X = {xi, • • ■ ,x„} , 
the convex hull Co(X) is defined as the minimal set 
containing all points in X , satisfying that Co (X) = 
{J2i=l a i X i \ x i £ X, a l > 0, Ya=i a i = 1 }■ 

Graph theory (cf. [26] and [27]) is widely used to rep- 
resent a networked system. Let Q = (V, £ ) denote a di- 
rected graph, where V = {v\, • • ■ , vn} and £ C V x V 
denote the set of nodes and the set of edges, respectively. 
Each edge (vi,Vj) G £ represents the neighborhood of 
node i and node j, which indicates that node i is able to 
access states of node j, but not vice versa. The neigh- 
bor set of node i is denoted as Afi = {vj | (vi,Vj) G £} . 
A directed path from node v\ to node Vk is a sequence 
of edges (i>i,i>2), {v%, V3) , • • ■ , (vi,Vk) in the directed 
graph Q. If graph Q contains a directed tree, every node 
has exactly one parent node except for one node, called 
the root, and the root has directed paths to every other 
node in graph Q. The adjacency matrix is defined as 
A = [a.j] e R NxN with a„ > 0 if {v l ,v J ) G £, and 
fly = 0 otherwise, where a,ij represents a weighting fac- 
tor for the associated edge (vi,Vj). A matrix with posi- 
tive or zero off-diagonal elements is referred to as a Met- 
zler matrix [28] . The Metzler matrix L for the graph Q is 



3 A continuous function a : [0, a) — > [0, 00) is said to belong 
to class K, if it is strictly increasing and a (0) = 0. It is said to 
belong to class ICoo if a = 00 and a (r) — > 00 as r — > 00 [22]. 



defined as L = A-D G R N x N , where D = [d^] G R NxN 
is a diagonal matrix with da = X^Li a ij- To facilitate 
the following development, a Lemma in [29] is introduced 
as follows. 

Lemma 2 [29] Consider a linear system x (t) = 
A (t) x it) , and a Lyapunov function V (x) = 

max j xi, • • ■ , x„ j — min j xi, • • • , x n } , where x (t) = 

[xi, . . . , x n ] T Gi? n is a n dimensional state. If the time- 
varying matrix A(t) G R nxn is a piecewise continuous 
function of time with bounded elements, A (t) is a Metzler 
matrix with zero row sums, and the time-varying graph 
corresponding to A (t) has a spanning tree for all t > 0, 
then V < 0 for allt > 0 and consensus is achieved, i.e., 
X\ • x n . 

3 Problem Formulation 

3. 1 Individual Social Behavior 

Consider a social network composed of n individuals. 
Each individual i maintains a state qt (t) G M d in a social 
network, which can represent opinions on social events, 
or human emotional states such as happiness, love, anger 
or fear. It is assumed that the current state qiit) of an 
individual i can be detected from other social neighbors 
such as close friends or family in the social network. Gen- 
erally, the opinions or emotional states formed by indi- 
viduals about social events are not only influenced by 
the information gathered through communication with 
their social neighbors, but also depend on their personal 
experiences. To capture the evolution of individual so- 
cial states by taking into account not only exogenous 
influence (e.g., information from friends or family) but 
also their own character (e.g., past experience, memory), 
inspired by the works of [8-11], qiit) is modeled as the 
solution to a fractional-order dynamics as 

0 D? qi (t) = Ui (t), i=l,.-.,n, (2) 

where m G K d denotes an influence (i.e., control input) 
over the social state, and oDfqiit) is the a th derivative 
of qi(t) with a G (0, 1]. 

Note that the model in (2) is a heuristic approximation 
to a social response, which indicates that a person's so- 
cial state has a direct relationship with external influ- 
ence integrated over the history of a person's emotional 
states. On-going efforts by the scientific community are 
focused on the development of clinically derived mod- 
els; yet, to date, there is no widely accepted model to 
describe a person's social response in a social network. 

3.2 Social Interaction 

Let Q (t) = (V, £ (t)) denote a directed graph, where the 
node set V = {!,■■■ ,n} represents individuals and the 
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edge set £ C V x V represents the interactions between 
individuals in a social network. Suppose that there exist 
to followers i G Vf, i — 1, • ■ ■ ,to, and n — to leaders 
j 6 Vl, j = m + 1, • • • , n, where Vl and Vf denote 
the leader and the follower set, respectively, satisfying 
Vl U Vf = V. It is assumed that the leaders' states 
q% G K d , i G Vl, are desired and immutable. For each 
follower, its state i G Vf, evolves according to the 
dynamics (2) under the influence from both followers and 
leaders directly or indirectly by the underlying network. 

A directed edge G £ in Q represents the neighbor- 
hood of node i and j. Each edge (i, j) is associated with 
a time-varying weighting factor called the social differ- 



ence Sij £ R , which is defined as SV 



9j\ 



Since individuals are assumed that they fail to incorpo- 
rate the information provided by neighbors whose states 
are far from their own, the designed social difference Sm- 
arms to capture the closeness of the states between two 
neighboring nodes i and j. It is also assumed that there 
exists a threshold S £ K + , and individuals i and j are 
able to influence each other only when their social dif- 
ference Sij < S. In other words, an edge (i, j) in graph Q 
does not exist if the social difference Sij is greater than 
the threshold 6. The neighbors of individual i in graph 
Q are defined as M% = {j \ Sij < 6}, which determines a 
set of individuals who can influence the social states of 
individual i. A directed path from node v\ to node Vk is 
a sequence of edges (v%, U2) , (V2, V3) , • ■ ■ , (vi, Vu) in the 
directed graph. If graph Q contains a directed spanning 
tree, every node has exactly one parent node except for 
one node, called the root, and the root has directed paths 
to every other node in graph Q . 

Assumption 1: For each follower i G Vf, there exists at 
least one leader that has a directed path to the follower 
i in the initial graph Q (0). 

Assumption 1 implies that there exists a directed span- 
ning tree for the initial graph G(0), where the set of 
leaders acting as the roots in the directed spanning tree 
has an influence directly or indirectly on all followers 
through a series of directed paths in the network. Note 
that a connected graph (i.e., a directed tree structure) is 
only assumed in the initial graph, and the controller de- 
veloped in the subsequent section will preserve the net- 
work connectivity to ensure consistent influence between 
social neighbors. 

3.3 Objectives 

Let q(t) , q F (<), and q L (t) denote the stacked vector 
of all states qi (t) , i G V, the followers' states qi (t) , 
i G Vf, and the leaders' states qi (t), i G Vl, respec- 
tively. The convex hull spanned by the states of lead- 
ers, and all states (i.e., both leaders and followers), are 
then represented as Co (q L (t)) and Co (q (t)) , respec- 
tively. Since the leaders' states are static, the convex hull 



Co (q L (£)) is constant, while the convex hull Co (q (t)) 
is time varying and depends on the states of the follow- 
ers. After formulating the social network as a networked 
fractional-order system described by (2), the objective 
is to regulate the states of followers to a desired region, 
which is a convex hull spanned by all stationary leaders' 
states (i.e., <&(£) — > Co(Vl) Vi G Vf). To ensure that 
each individual is able to be influenced by social leaders 
through a path of directed edges by communication with 
their local neighbors only, another goal is to preserve the 
network connectivity for the underlying social network 
(i.e., maintain the social difference Sij < 6 so that peers 
remain peers) when given an initially connected graph 
Q. Since the systems in (2) along different dimensions 
are decoupled, for the simplicity of presentation, only a 
scalar system (d = I), that is qi(t) G K, is considered 
in the following analysis. However, the results are valid 
for a d dimensional case by the introduction of the Kro- 
necker product. 

4 Distributed Influence Design 

The artificial potential field based approach is one of the 
most widely used methods in the control of multi-agent 
systems, which consists of an attractive potential encod- 
ing the control objective and a repulsive potential rep- 
resenting the motion constraints (cf. [30]). To apply the 
potential field based approach to a social network prob- 
lem, inspired by the work of [31] and [32], a decentralized 
potential function ipi : R d — > [0, 1] Mi G Vf is developed 
to influence the followers' states to a desired end as 



7i 



7f + 



A) 



i/fe ' 



i G Vf 



(3) 



where k G 



is a tuning parameter, ji 



is the goal function, and (3i : W 1 — > 1R + is a constraint 
function. The goal function in (3) is designed as 



H = 



(4) 



which aims to achieve consensus for node i with its neigh- 
bors j G Mi . To ensure consistent influence from neigh- 
bors (i.e., maintain the social difference Sij < 5), the 
constraint function in (3) is designed as 



Pi 



where b. 



ij 



5 - Sa G 



(5) 



". For an existing interaction 



between individuals i and j, the potential function ipi in 
(3) will approach its maximum whenever the constraint 
function /3j decreases to 0 (i.e., the social difference S^ 
increases to the threshold of 8). 

Based on the definition of the potential function in (3), 
the distributed influence algorithm for each follower is 
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designed as 

u t = -KiS7 qi (pi, i G V F 



(6) 



where Ki is a positive gain, and V 9i tfi denotes the gra- 
dient of (fi with respect to c^. Applying (6) to (2), the 
closed-loop dynamics of social response for all individu- 
als in a social network can be obtained as 



oD? qi (t) = -KiV qi (pi ieV F 
0 D? qi {t)=Q iEV L . 



(7) 



Since leaders' states are stationary, the input to leaders 
in (7) is zero, and S7 qi ifi can be computed as 



V 9i ipi 



(8) 



From (4) and (5), V Qi ji = J2jeJVt (<H _ % 0 and 
= - SieM &y (© _ Sj) respectively, where by = 
ILeA/i.J^i ^ e R+ ' Substituting V 9i 7i and V 9i /3j into 
(8), V qi (pi is rewritten as 



V 9i </>i = X] - eM (® - «i) j 



where 



fc(7? +A) 



(9) 



(10) 



is non- negative, based on the definition of 7$, /3j, fc, and 



5 Convergence Analysis 



do not have a directed path to influence the follower's 
states. Hence, the state of the follower is independent 
of the influence of leaders, and thus can not converge to 
the stationary convex hull spanned by leaders. To en- 
sure the continued influence from leaders to all followers, 
a directed spanning tree structure must be maintained 
all the time. The following theorem shows that, given 
an initial graph containing a directed spanning tree as- 
sumed in Assumption 1, the tree structure will be pre- 
served under the influence function in (6) (i.e., network 
connectivity is maintained and social peers do not be- 
come isolated from the social group). 

Theorem 1 The influence function in (6) guarantees a 
directed spanning tree structure in Q for all time. 



PROOF. It is assumed that the initial graph Q (0) has 
a directed spanning tree in Assumption 1. If every ex- 
isting edge in Q (0) is preserved, the tree structure will 
also be preserved. Since an individual only considers a 
local neighbor's opinion as reasonable when their social 
difference Sij < 5, peer influence is maintained when 
each edge Sy < 6 all the time. Consider a state qo for 
individual i, where the interaction between individual i 
and neighbor j G A/i satisfies 6y (go, <y) = 0, which indi- 
cates that their social difference is too large to influence 
each others' opinion, and the associated edge is about to 
break. From (5), /3; = 0 when 6y = 0, and the naviga- 
tion function ifi achieves its maximum value from (3). 
Since ipi is maximized at qo , no open set of initial con- 
ditions can be attracted to qo under the negative gra- 
dient control law designed in (6). Therefore, the social 
bond between individual i and j is maintained less than 
5 by (6), and the associated edge is also maintained. Re- 
peating this argument for all pairs, every edge in Q is 
maintained and the directed spanning tree structure is 
preserved. □ 



To show that the followers in the fractional-order net- 
work converge to a convex hull spanned by the static 
leaders' states, the following analysis is segregated into 
three proofs. The first proof shows that the existing in- 
teraction between individuals is maintained by the influ- 
ence function designed in (6) (i.e., the social difference 
Sij < 5 for all time, meaning influential peers remain 
influential), and thus the connectivity for the social net- 
work is preserved. The second proof yields the asymp- 
totic stability for an integer-order representation of the 
dynamic system in (2), which is then used to establish 
the asymptotic convergence to the equilibrium set of con- 
sensus states for the fractional-order system by using a 
Mittag-Leffler stability analysis in the third proof. 

5. 1 Maintenance of Social Influence 



5.2 Convergence Analysis 

To establish asymptotic convergence to the equilibrium 
points (i.e., the convex hull Co (q L (£)) for the fractional- 
order dynamics in (2), an integer-order system qi(t) = 
Ui(t) with a = 1 in (7) is considered first in the following 
theorem. 

Theorem 2 Consider a network composed of station- 
ary leaders and dynamic followers described by (7). The 
followers i G Vf asymptotically converge to the equilib- 
rium points (i.e., a convex hull Co (q L (t)j spanned by 
the leaders states only), if there always exists at least one 
leader j G Vl that has a directed path to any follower i 
(i.e., a directed spanning tree is maintained). 



If a directed graph Q (t) does not have a directed span- 
ning tree, there must exist a follower to which all leaders 



PROOF. This theorem is proven by using LaSalle's in- 
variance principle and convex properties. Let V(q(i)) 
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be the volume of the convex hull Co (q (t)) formed by all 
leaders' and followers' states. First, we show that there 
exists a compact set f2 such that if qi (0) G f2 for Vi G V, 
then qi (t) G il for all £ > 0, which implies that f2 is a 
positively invariant set. Second, let E be the set of all 
points in f2 where V = 0 (i.e., the volume of Co (q (£)) , 
q(£) G E, stays constant). It is then shown that M is 
the largest invariant set, where M is the set of points in 
the convex hull Co (q L (£)) formed by stationary leaders 
only. 

Substituting (9) into (7) with a = 1 yields the following 
closed-loop emotion dynamics 

= - 2 K i m tj (t) - ?j (*)) i g Vf 

j'GA/i (11) 

&(f) = ° J6V L , 

which can be rewritten in a compact form of a time- 
varying linear system as 



q(*) 



where 0 







Tr(i) 


A L (*)_ 




_ ^(n— m) x n _ 



q(f)> 



(12) 



(n-m)xn denotes the (n — m) x n matrix with all 

mx ™ are defined as 



zeros, and the elements of 7r (£) G K 

i = k 



' i j 



- Ki m 

n ik (t)={ K . m . k k GAfi,i^k 



0, 



k£Ni,i^ k. 



Each follower i 6 Vf in (11) evolves according to the 
dynamics: 



(14) 



To facilitate the analysis, (14) can be written in discrete 
time as 




The next step is to show that all followers' states will 
asymptotically converge to their equilibrium points. To 
see that the equilibrium points are indeed the station- 
ary convex hull Co (q L (£)), let qi_ eq be an equilibrium 
point for a follower i 6 Vf- For an equilibrium point, 
it must have qi >e q = 0, and (14) can be written as 
_ Y^jeMi n iJ (*) (9%,eq (*) - Qj,eq (t)) = 0, which yields 
that 



1 



(16) 



by using (14). Due to the fact that —ttu (£) = 
EjeA/< ( ) ^ ^ + n " om (13), (16) indicates that the 
equilibrium point g,, eg lies in a convex hull spanned by its 
neighbors' states (i.e., leaders and/or followers). Since 
every follower ends up in a convex hull spanned by its 
neighbors' states and the leaders' states are static, every 
follower will end up in a convex hull spanned by the lead- 
ers states only (i.e., Co (q L (£))). Using the fact that 
is non-negative from (10) and K{ is a positive constant 

TT(t) 



gain in (6) , 



0 



(n—m) xn 



in (12) is a Metzler matrix with 



zero row sums. According to Lemma 2 and following a 
similar procedure as in [29], the convex hull Co (q (£)) is 
shrinking (i.e., V (q (£)) < 0), since the difference of the 

extremes max (xi, • • • , x„ j and min | x\, • • • , x n j 
is decreasing. If all followers states are initially within 

the states will always stay 
0). 



the convex hull Co (q L (£)) 
(13) within Co (q L (£)) (i.e., V 



A proof by contradiction can now be used to show that M 
(i.e., Co (q L (£))) is the largest invariant set. Let M' D 
M be a larger invariant set in E. Suppose that there is 
a follower whose state qi (0) ^ M, and qi (0) is on the 
boundary of M' . Since M' C E, the volume of the set 
M' stays constant. The only way for the volume of M' to 
stay constant is that qi (0) = qi (t) for all t > 0. However, 
for this to happen, we must have mj (£) = 0 for Vj G A/i 
from (14), which indicates that the follower i is isolated 
from the group. This isolation is a contradiction with 
network connectivity. Hence, M is the largest invariant 
set. The followers asymptotically converge to the largest 
invariant set M (i.e., the equilibrium points Co (q L (£))) 
T^ijQj (t) , (15) by using LaSalle's invariance principle in [22] . □ 



where T is a sufficiently small sampling period. From 
(15), it is clear that <&(£ + 1) is a convex combination 
of its current state qi (t) and its neighbors' states qj (t), 
j G Afi, which implies that the follower i moves to- 
wards the convex hull spanned by itself and its neighbor- 
hood set Afi. Since the leaders' states are stationary and 
the followers' states are evolving within the convex hull, 
V (q (£)) is uniformly non-increasing and thus V (q (0)) 
is the compact set O. 



Since asymptotic stability for the integer-order system 
(11) is established in Theorem 2, a similar proof proce- 
dure in our recent work [18] can be followed to prove 
asymptotic stability for the fractional order system in 
(7) by using Mittag-Lefher stability analysis and a Con- 
verse Lyapunov Theorem. 

Theorem 3 The follower i G Vf, with closed-loop 
fractional- order dynamics in (7) with a G (0, 1), asymp- 



6 



totically converges to the convex hull spanned by station- 
ary leaders if at least one leader j G Vl has a directed 
path to the follower i. 

PROOF. Let Xl (t) = q, (t) - q i>eg , and x F (t) = 
q F (t) — q F , where q F denotes the stacked vector of 
q% t eq- Since the leaders' states are constant, the closed- 
loop fractional-order dynamics in (7) can be written in 
a compact form as 

0 D^ F (t) = , 9 (x F ), (17) 

for all followers where g(x F ) is a function of follower 
states. Since stability of a fractional-order system is de- 
fined by Definition 1, and Mittag-Leffler stability implies 
asymptotic convergence as discussed in [24], the follow- 
ing development aims to show that (17) is Mittag-Leffler 
stable. 

Since asymptotic stability is established in Theorem 2 
for the integer-order system of (15), a Converse Lya- 
punov Theorem, (i.e., Theorem 4.16 in [22]) is invoked 
to establish that there exists a function l/(t,x F ) : 
(0,oo] x R m ^ R and class JC functions a t (i = 1,2,3) 
such that 

^i(||x F ||) <V(t,x F ) <CT 2 (||x F ||), (18) 
V < -OS (||x F ||) . (19) 

Let /3 = 1-a G (0, 1). From Theorem 8 in [24] and (19), 
the fractional derivative of V is computed as 

0 D?V {t, x F ) = Q D\-°V (t, x F ) = 0 D- a V (20) 
< - 0J D i - a (a 3 (||x F ||)). (21) 

From the definition of the fractional integral 
oA" Q /W = rfey/o (Jff-^ T, where T (•) denotes 
the Gamma function [7], it is known that oD^~ a (■) is a 
class JC function, since oD^ a (0) = 0 and oD^ a (■) is 
strictly increasing on the domain [0, oo). Using the fact 
that ai o a 2 also belongs to class /C, where a\ and 02 are 
class JC functions, the term oD^ a (03 (||x F ||)) in (21) 
is a class JC function, since both oD^ a (•) and 173 (•) are 
class JC functions. Thus, the inequality in (21) can be 
written as 

0 D?V(t,x F ) <-a 4 (||x F ||), (22) 

where 04 (||x F ||) = oDj~ a (0-3 (||x F ||)) is a class JC func- 
tion. Applying Lemma 1 to (18) and (22), Mittag-Leffler 
stability of (17) with a G (0, 1) can be obtained, which 
implies that the equilibrium points q F for the follow- 
ers in the closed-loop fractional-order system in (17) are 
asymptotically stable. □ 
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Figure 1. The Zachary's karate club network is modeled by 
a directed graph. 

6 Simulations 

To illustrate the proposed control algorithms, simula- 
tions are performed on a karate club network described 
in [33]. The karate club network considered in this sec- 
tion consists of 3 social leaders and 7 followers, and is 
represented as a directed graph in Fig. 1. The solid ar- 
row connecting two individuals in Fig. 1 indicates an 
established social bond (e.g., friendship) and the di- 
rected influence between individuals. Note that the lead- 
ers can not be influenced, while the followers can be in- 
fluenced by social peers as well as social leaders. Based 
on the topology described in Fig. 1, each individual is 
randomly assigned a social state (e.g., an opinion on 
an event or an emotional state such as happiness, fear 
and anger). Without loss of generality, we assume that 
the social states of individuals are two dimensional (i.e., 
qi(t) Gl 2 ). The control law in (6) yields the simulation 
results shown in Fig. 2, which illustrates that the fol- 
lowers' states converge to the convex hull formed by the 
social leaderfH 

7 Conclusion 

By modeling the group social response as a networked 
fractional-order system, a decentralized potential field- 



4 An infinite state model is developed in [25] to represent a 
Fractional-order Differential Equation (FDE), as a means to 
solve the initial condition challenges associated with FDEs. 
As indicated in [25], the internal state z (t) in the infinite 
state model contains historical information of the fractional- 
order system. To address the initialization problem of a FDE, 
the observer based model in [25] can be used to estimate 
z (to ) by using past information from t < to . Since the current 
work focuses on showing how the individual's behavior can 
be influenced under the designed influence function, for sim- 
plicity in the simulation, z (to) is assumed known, which cap- 
tures the individual's historical experience. Given that the 
fractional-order dynamics in (2) is written as qi (t) = I n (m) 
where (114) is the nth fractional integral of control input Mi, 
the trajectory of qi (t) is simulated by following the infinite 
state approach in [25]. An alternate initialization approach 
is to account for the initialization function (cf. [34,35]). 
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Figure 2. Evolution of social states. The followers' states 
converge to the convex hull formed by social leaders, where 
the the leaders' social states are denoted as squares, and the 
followers' initial and final social states are denoted by circles 
and dots, respectively. 

based influence algorithm is developed in this work to en- 
sure that all individuals' states achieve consensus asymp- 
totically to a desired convex hull spanned by the sta- 
tionary leaders' states, while maintaining consistent in- 
fluence between individuals (i.e., network connectivity). 
This work considers individuals whose social response is 
modeled by a FOS with a € (0,1]. Since some individuals 
may respond with a more complex dynamic (e.g., a 6 (1, 
2]), future efforts will focus on generalizing the develop- 
ment to include networks with heterogeneous members 
with higher order dynamic response. Future effort will 
also consider different influence capabilities between in- 
dividuals. For instance, a person tends to have a larger 
tolerance for a difference of opinions for a certain social 
event in a close friend than a loose acquaintance, and 
thus, can be more easily influenced by the close friend. 
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